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ABSTRACT 

We examine the chemical and dynamical structure in the solar neighbourhood of a 
model Galaxy that is the endpoint of a simulation of the chemical evolution of the 
Milky Way in the presence of radial mixing of stars and gas. Although the simulation's 
star-formation rate declines monotonically from its unique peak and no merger or 
tidal event ever takes place, the model replicates all known properties of a thick disc, 
as well as matching special features of the local stellar population such as a metal- 
poor extension of the thin disc that has high rotational velocity. We divide the disc 
by chemistry and relate this dissection to observationally more convenient kinematic 
selection criteria. We conclude that the observed chemistry of the Galactic disc does 
not provide convincing evidence for a violent origin of the thick disc, as has been 
widely claimed. 

Key words: galaxies: abundances - galaxies: evolution - galaxies: ISM - galaxies: 
kinematics and dynamics - Galaxy: disc - solar neighbourhood 



1 INTRODUCTION 

Our Galaxy's stellar disc was first divided into two compo- 
nents because the vertical density profile derived from star 
counts could be fitted by a superp osition of two exponen - 
tials but not by a single exponential (|Gilmore &; Reidlll983T ). 
Further investigations revealed a thick-disc component that 
was characterised by a high velocity dispersion, high a en- 
richment and a remarkably old age. Many authors consider 
the thick disc to be a relic of a turbulent era of Galactic 
history in which the thick disc formed from accreted satel- 
lites and/or a thin disc was violently heated by one or m ore 
merger events (for a discussion see e.g. lReddv et al.ll2006h . A 
violent origin of the thick disc is strongly supp orted by tradi- 
tiona l models of chemical evolution (see e.g. lChiappini et all 
I1997T ). These models require a period of rapid star forma- 
tion early in the life of the disc, followed by a period in 
which star formation effectively ceased in which the inter- 
stellar medium could be enriched with iron by SNIa and the 
overall metallicity level could be brought down by accretion 
of metal-poor gas. 

In an earlier paper (Schonrich & Binney] 2009, here- 
after SB09) we showed that when one allows for radial mix- 
ing, which is an unavoidable c onsequence of spir al struc- 
ture (|Sellwood fe Binnevlliooi iRoskar et al.ll2008l ). a two- 
component disc arises naturally in the simplest model, in 
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which the star- format ion rate (SFR) is a monotonically de- 
clining function of time from its global maximum. In this pa- 
per we examine in greater depth the contents of the model's 
solar neighbourhood, especially its characteristic stellar pop- 
ulations. We identify the solar-neighbourhood signatures of 
the thin and thick discs and analyse the relationship be- 
tween the chemistry and kinematics of nearby stars. This 
exercise enables one to understand better the relationship 
between the underlying nature of the thin and thick discs 
and samples of stars that have been selected by kinematic, 
chemical or spatial criteria. A better understanding of this 
relationship is of considerable practical importance because 
substantial allocations of telescope time are currently com- 
mitted to spectroscopic surveys (SEGUE, RAVE, HERMES, 
WFMOS, Gaia) that are designed to unravel the nature and 
history of the thick disc, and a clear picture will not emerge 
from these surveys without a secure understanding of selec- 
tion effects. 

The paper is structured as follows. In Section 2 we sum- 
marise the physics that underlies the model and analyse its 
prediction for the disposition of solar-neighbourhood stars 
in the ([Fe/H], [a/Fe]) plane. We identify the thin and thick 
discs within this plane, and show how the kinematics and 
ages of stars vary within the ([Fe/H], [a/Fe]) plane. The 
structures we identify seem to have all the properties ex- 
pected of the Galaxy's thin and thick discs. In Section 3 we 
explore the extent to which stars can be successfully assigned 
to the thin and thick discs by kinematic selection. We show 
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Figure 1. The distribution of birth radii of stars in the model 
GCS sample (green dashed line) and of all stars in the solar cylin- 
der (solid red line). 



that this process cannot be clean. Finally in Section 4 we 
sum up and relate the characteristics of the features we have 
identified to the formation history of our model. Since the 
SFR in the model has been monotonically decreasing from 
its global maximum, and no merger or tidal event has ever 
occurred in it, we conclude that, contrary to widespread be- 
lief, features of the Galaxy, such as the overlap in [Fe/H] of 
the thin and thick discs do not in fact constitute convincing 
evidence for a violent origin of the thick disc. 



2 THE MODEL 
2.1 Physical inputs 

The SB09 model is the endpoint of a simulation of chemi- 
cal evolution within a disc in which the star- format ion rate, 
whic h is controlled b y the surface density of the ISM as in 
the lKennicuttl (| 19981 ) law, declines monotonically with time 
from a unique global maximum. The gas disc always has 
an exponential surface density with scale length 3.5 kpc so 
that by the Kennicutt law the young stellar disc has an 
exponential surface density with scale length 2.5 kpc. The 
assumptions regarding the (universal) initial mass function, 
stellar lifetimes and yields are also taken from the litera- 
ture. The only novel features are a radial flow of gas within 
the disc and radial migration of stars. The latter occurs 
both because over time stars change their angular momenta 
( "churning" ) and because they move on orbits that become 
increasingly eccentric and inclined to the Galactic plane 
("blurring"). Traditional models of Galactic evolution have 
ignored these effects although it has always been evident 
that blurring occurs. The importa n ce of churning was only 
realised when ISellwood &; Binnevl (2002) found that even 
weak spiral structure in N-body simulations causes stars to 
shift their guiding centres by a kiloparsec and more in a 
single dynamical time. These motions, which arise when a 
star is at the corotation resonance with a spiral arm, do not 
heat the disc, so they come to light only when the angu- 
lar momenta of individual stars are followed. In the model 
the extent of churning is governed by a parameter k c h that 
could be determined from N-body models if we knew the 
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Figure 2. The metallicity distribution of solar-neighbourhood 
stars: data points from Holmberg et al. (2007); red curve the SB09 
model. For the model total metal abundance is plotted horizon- 
tally, while for the data the plotted quantity is the photometric 
metallicity indicator given in Holmberg et al. (2007). 



past strength of spiral structure. SB09 determined k c h by 
fitting the model to the distribution of solar-neighbourhood 
stars in [Fe/H]. The radial dependence of churning strength 
was taken to be proportional to the product of surface den- 
sity and radius Si?, following an argument based on disc 
instabilities. 

The dashed green line in Fig. [T] shows the distribution 
of birth radii of stars in the model GCS sample. Because 
the GCS stars lie near the plane, the fraction of these stars 
that are young is higher than the fraction of young stars 
in the entire solar cylinder. This bias towards young stars 
leads to the distribution of birth radii of GCS stars being 
narrower than the distribution of birth radii for all stars 
in the solar cylinder, which is show by the full red curve 
in Fig. [U The difference between the two distributions is 
largest for stars born at < 5 kpc because those inner disc 
stars have larger vertical velocity dispersions and therefore 
larger scaleheights. 

Hot gas was assumed to be too far from the disc to take 
part in churning, while the cold gas and stars were assumed 
to be equally involved in this process. It is likely that these 
assumptions exaggerate the impact of churning on old stellar 
populations, which have high velocity dispersions, relative to 
its impact on young stellar populations. Since there is as yet 
no basis for quantifying the impact of velocity dispersion 
on churning rate, the model of Paper I avoids additional 
undetermined parameters by ignoring this possibility. 

The flow of gas within the disc enables the surface den- 
sity of star formation to be an exponential function of ra- 
dius even though the rate at which accreted gas joins the 
disc does not necessarily vary exponentially with radius. The 
surface density of inflow of metal-poor gas and the flow of 
gas within the disc are jointly controlled by two parame- 
ters, /a and /b, which substitute for a knowledge of the 
radial profile of cosmic infall. Attempts to o btain the latter 
from simulations (e.g., IColavitti et alll2008l ) have not been 
successful, probably because much of the gas that joins the 
disc spends a significant time after infall in the warm-hot in- 
tergalactic medium (WHIM). Hence at the present time we 
must parametrise the infall in some way. SB09 found that 
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/b is effectively set by the measured oxygen gradient in the 
ISM, and /a was fitted alongside k c h to the local metallic- 
ity distribution of stars. Fig. [2] shows the fit that was ob- 
tained to data for ~ 10 000 non-binary stars in the Geneva- 
Cope nhagen survey (|Nordstrom et al. 120041 ; iHolmberg et al.l 
l2007l . hereafter GCS). 

The random velocities of stars formed at a given 
radius are assumed to increase with age r as r 1 ^ 3 in 
line with the predictions of both theory (Ijenkinsl 1992) 
and st udies of the sola r neighbourhood (|just &; Jahreissl 
20071: lAumer fe Binnevl 120081 ). We have modified the 
Schonrich &; Binnevl (|2009) model very slightly by increas- 



ing the (vr) 1 ^ of a 10 Gyr-old popula tion of local stars from 
38kms _1 fcf lDehnen fe Binnevlll998 ) to 45kms~ 1 . This in- 
crease brings the (vr) 1 ^ 2 for the entire GCS sample into 
line with the observed value. In c oncordance with observa- 
tions (e.g. iLewis &; Freeman! 1 19891 ) the square of the intrin- 
sic velocity dispersion at a given age is assumed to scale 
with radius as e -- R / 1 - 5 - R * ) where R* = 2.5 kpc is the scale- 
length of the stellar disc, so (vr) 1 ^ 2 ~ 90kms _1 at R = R*. 
The square of the vertical velocity dispersion component is 
assumed to show a slightly steeper rise (implying approxi- 
mately constant scaleheight) being proportional to e ~ R / R * . 



2.2 Predictions of the model 

The merit of the SB09 model is that it tracks the kinematics 
of stars in addition to their chemistry. Observations always 
have a kinematic bias of some kind, either because they are 
restricted to stars that lie near the Sun and therefore the 
plane, a region favou red by stars with small vertical velo city 
dispersions, or (as in ljuric et alfeooil : llvezic et al.ll2008l ) be- 
cause they focus on faint stars that are far from the plane, or 
because an explicit high- velocity criterion is applied in order 
to reduce contamination of a thick-disc sample by thin-disc 
stars. A model that includes both chemistry and kinematics 
is essential for the interpretation of a kinematically selected 
study. 

The SB09 model makes predictions for the global struc- 
ture of the Galaxy's stellar and gas discs, but in this pa- 
per we focus on the solar neighbourhood, and especially 
the stars that happen to lie within 100 pc of the Sun. This 
volume is of particular interest because within it G-dwarfs 
are bright enough to have good Hipparcos parallaxes, mea- 
sured radial velocities, and in a few cases medium-to-high 
resolution spectra from which detailed chemical abundance 
patterns can be extracted. The GCS provides space ve- 
locities, surface gravities and metallicities for this volume, 
and detailed abundance analyse s have been carried out for 
much smaller subsets o f stars (lFuhrmannll998l;lBensbv et al.l 
2003 ;IVenn et al.ll2004l : iBensbv et al.ll2005l : iGilli et alJfeOOd 
Red dv et al.l [2006) . Because the GCS sample is essentially 
magnitude limited, it is not representative of the volume 
typically simulated by models of chemical evolution, namely 
a cylindrical annulus around the disc. In particular, thick- 
disc stars are under-represented within the GCS relative to 
a cylindrical annulus. The SB09 model provides an arena in 
which the impact of this bias can be assessed. 

Fig. [3] shows the densities of simulated GCS stars in 
the ([Fe/H], [Ca/Fe]) (upper) and ([Fe/H], [O/Fe]) planes. 
Trajectories of the cold ISM at galactocentric radii of 
10, 7.5, 5, 2.5 kpc (from left to right) are indicated by black 




Figure 3. Logarithmic stellar densities for a simulated GCS stel- 
lar sample in the ([Ca/Fe], [Fe/H]) (top) and ([O/Fe], [Fe/H]) 
planes. Contours have a 0.5 dex spacing. Black lines track the 
development of the cold ISM in annuli of radii of (from right to 
left) 2.5,5.0,7.5 and 10 kpc. 



lines. The ISM starts at early times with low metallicity and 
high a enhancement at the top left of each panel. As the gas 
is enriched with metals, each trajectory moves to the right. 
With the onset of SNIa, the composition of stellar yields 
shifts towards iron, so [a/Fe] decreases and the trajectories 
move downwards. Eventually the ISM approaches a steady 
state in which additional enrichment is balanced by the infall 
of fresh metal-poor material from the IGM. Since the delay 
of SNIa-enrichment is assumed to be independent of envi- 
ronment, the time at which trajectories first move downward 
is independent of radius. By the fact that the timescale of 
SNIa-enrichment does not vary with radius either, the ISM 
trajectories tend to be nearly aligned. Thus the point of 
turndown is at higher metallicities (further to the right) for 
populations closer to the Galactic centre, where the ISM is 
enriched faster by more intense star formation relative to 
the present gas mass. 

The colours and green contours in Fig. [3] show the den- 
sity of stars within each plane. In each panel two ridges 
of high density are apparent - one at high [a/Fe], which 
we call the metal-poor thick disc and one at low [a/Fe], 
which is associated with the thin disc. Crucially, the thin- 
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disc ridge runs at a large angle to the black trajectories of the 
ISM. Thus the thin-disc ridge in no sense traces the chemi- 
cal history of the thin disc; instead it reflects the spread in 
radii of birth of local stars, which gives rise to a spread in 
[Fe/H] by virtue of the metallicity gradient within the ISM 
(which is larger in the SB09 model than in traditional mod- 
els of chemical evolution). In a similar manner the thick-disk 
ridge follows the evolution of all rings at low metallicities, 
but stretches significantly to higher [Fe/H] than the point 
at which the solar annulus leaves it. 

The depression in the stellar density between the ridges 
in Fig. [3] is a consequence of the rapid downward motion 
of all trajectories after the onset of SNIa, and of the ISM 
approaching a steady state as it enters the thin-disc ridge 
line; relatively few stars are formed at intermediate values 
of [a/Fe]. Variations in the timescales of enrichment will 
change the depth of the depression - for example, a shorter 
timescale for the decay of SNIa progenitors will cause tra- 
jectories to move downwards faster, leading to fewer stars 
in the intermediate region. Our models use a prescription 
for SNIa that is standard for models of chemical evolution, 
with no SNIa until 0.15 Gyr after star formation, and then 
an expo nential decay in the rat e of SNIa with time constant 
1.5 Gvr. iMannucci et al.1 (|20Q6h suggest that -half SNIa ex- 
plode promptly (within 0.1 Gyr of star formation) and the 
rest explode at a rate that declines exponentially with time 
constant 3 Gyr. The existence of prompt SNIa would not 
materially affect our work as long as a significant fraction 
of SNIa are in the population with a long time constant, 
since the prompt SNIa will lower the alpha-enhancement 
level of the thick disc component, but not affect the evolu- 
tion between t he tw o density ridges in the abundance plane. 
iForster et all (2006) showed that timescales are very weakly 
constrained by SNIa counts due to uncertainties in the star- 
formation histories. 

Since both Ca and O are a-elements, the distributions in 
the upper and lower panels of Fig.[3]are qualitatively similar, 

and O will be the only a-element explicitly discussed below. 

I 1 \ — — \\ 

As lHavwoodl (2008) has pointed out, the principal trac- 
ers of radial mixing are the large dispersion in the metal- 
licities of stars in the solar neighbourhood and the strong 
increase in this dispersion with age, which is caused by im- 
migration of high-metallicity stars from inwards and low- 
metallicity stars from outwards. In fact, as SB09 demon- 
strated, it is impossible to fit the shape of the local metallic- 
ity distribution under plausible assumptions without radial 
mixing. 

As well as generating a large dispersion in the metal- 
licities of old stars, radial migration has a big impact on 
the interdependence of kinematics and chemistry. This im- 
pact is illustrated by Fig. 2] which is another plot of the 
([Fe/H], [O/Fe]) plane, but now with colour indicating the 
mean rotation velocity of stars at each point - blue indi- 
cates low rotation velocities and red large ones. We see that 
at any given value of [a /Fe] , there is a tight correlation be- 
tween [Fe/H] and rotation velocity in the sense that high 
[Fe/H] implies low rotation velocity because stars with high 
[Fe/H] are migrants from small radii and tend to be defi- 
cient in angular momentum, while stars with low [Fe/H] are 
migrants from large radii. The black lines that show the tra- 
jectories of the ISM almost constitute contours of constant 
mean rotation velocity, but there is a barely perceptible ten- 
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Figure 4. The structure of a simulated sample of GCS stars in the 
([Fe/H], [O/Fe]) plane. Contours spaced by 1 dex give the density 
of stars in this plane, while colour codes the average rotational 
velocity of the stars found at the point in question - the local cir- 
cular speed is assumed to be 220kms — 1 . Black lines give the tra- 
jectories of the cold ISM during the model Galaxy's evolution for 
galactocentric distances of (from left to right) 10, 7.5, 5, 2.5 kpc. 



dency for the rotation velocity to decrease as one moves up 
along a black line. 

While the correlation between [Fe/H] and rotation ve- 
locity seen in Fig. [4] is qualitative consistent with stars being 
scattered to more eccentric orbits while retaining their angu- 
lar momenta ( "blurring" ) , quantitatively changes in angular 
momentum ( "churning") play a big role in structuring Fig. 2] 
While churning does occasionally move the guiding-centre 
radius of a star from R g < Ro to R g > Ro and thus increase 
the mean rotation speed at large [Fe/H], the dominant ef- 
fect of churning is to move guiding centres from R g <C Ro 
to R g < Ro such that a very metal-rich star is found in the 
solar neighbourhood at a relatively low rotation velocity. To 
illustrate the impact of churning quantitatively, if angular 
momentum were conserved, the population born 5 kpc from 
the Galactic centre would have ~ 150kms _1 , while stars 
born at 10 kpc would have ~ 300kms _1 . In the simu- 
lated sample, the mean speeds associated with these radii of 
birth are actually 200kms _1 and 240kms _1 . 

Churning substantially increases the chemical hetero- 
geneity of the stars that one finds near the Sun with a given 
V velocity: if high-a, high- [Fe/H] stars were brought to the 
Sun only by blurring, then all stars with a given V and there- 
for angular momentum would have identical chemistry. By 
changing the angular momenta of stars, churning ensures 
that stars of a given chemical composition are seen near the 
Sun over the whole range in V. 

In Fig. [4] the density of stars is indicated by white con- 
tours, which are spaced by 1 dex. The very top edge of the 
populated region is shaded blue, indicating low rotation ve- 
locities. Thus the highest-a stars form a structure with a 
large asymmetric drift. This fact reflects the large velocity 
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Figure 5. Same as Fig. [4] with colour coding for age and 0.5 dex 

spacing for density contours. 

dispersion of these stars, and the increased opportunity for 
migration enjoyed by this old population. Note that in this 
region of the diagram the colour rapidly changes to red as 
one moves downwards. Hence there is also a population of 
high- a stars that have large rotation velocities. These are 
typically slightly less old stars that formed outside the solar 
radius. The thin-disc ridge line at low a enhancement ranges 
from the strongly trailing at high metallicities to high rota- 
tional velocities and slig ht a enhan c emen t at low [Fe/H], as 
found observationally by iHavwoodl (|20Q8h . 

In Fig. [5]the contours show stellar density in increments 
of 0.5 dex and colours encode the mean age of stars, with 
blue implying youth. Naturally the oldest stars are high up, 
in the region of high [a/Fe]. Right at the top, lines of con- 
stant age run almost horizontally. As one descends the dia- 
gram, lines of constant age slope more steeply down to the 
right as a result of the more rapid decline in [a/Fe] at small 
radii. The youngest stars both from outer and from inner 
rings have yet to reach the solar neighbourhood in signifi- 
cant numbers, so in Fig. [5] several white contours are crossed 
as one moves along the thin-disc ridge from the location of 
the solar-radius ISM. 

In Fig.[6]both colours and contours (lOkms -1 spacing) 
show the velocity dispersion au of stars of a given chemical 
composition. The structure of the figure is the product of two 
mechanisms: (i) The velocity dispersion of stars born at any 
given radius scales with the third power of age, so older stars 
have larger velocity dispersions than younger stars born at 
the same locations. Consequently, in the figure velocity dis- 
persion tends to increase from bottom to top. (ii) Velocity 
dispersion increases inwards, so stars that have reached the 
solar neighbourhood from small radii of birth have larger 
velocity dispersions than stars that have reached us from 
large radii. When this fact is combined with the fact that 
for any given date of birth more metal-rich stars are born at 
smaller radii, a steep rise in au from left to right arises in 
Fig. [6] This plot suggests that we should include the region 
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Figure 6. Velocity dispersions (in kms -1 ) as functions of posi- 
tion in the ([Fe/H], [O/Fe]) plane. The graph is derived for a solar- 
neighbourhood sample by measuring the velocity dispersions of 
the populations with a specific chemical fingerprint. Two effects 
act on the velocity dispersion: The dependence on age mostly 
induces a top-down gradient following the evolution lines of the 
ISM. In the perpendicular direction (left-right) velocity dispersion 
increases with decreasing galactocentric radius. The low disper- 
sion of the Galactic thin disc is visible on the lower left side, 
girded by a high dispersion band running from top left to bottom 
right. 



of high velocity dispersion along the upper right part of the 
populated region in the thick disc. 

Since the single populations have - according to their 
ages and places of birth - different vertical dispersions, the 
older populations and those coming from inner radii will 
have higher scale heights and so reduced weights in a local 
sample. However, these populations dominate the compo- 
sition high above the Galactic plane. The upper panel of 
Fig. [3 depicts the iron abundance distributions of the stars 
at different heights above the plane. Both tails of the distri- 
bution are strengthened as one moves away from the plane. 
The growth in the proportion of metal-rich stars with \z\ 
is at first unexpected, but is a natural consequence of the 
higher vertical velocity dispersion of stars in the inner disc. 
Notwithstanding the growth of the metal-rich wing of the 
metallicity distribution, the mean metallicity falls with in- 
creasing \z\ by more than 0.2 dex, while the dispersion in- 
creases from below 0.3 dex to 0.5 dex. We expect the model, 
however, to underestimate the vertical metallicity gradient 
on account of our assumption that a star inherits the velocity 
dispersion of the galactocentric radius at which it was born. 
A better model would take account of the actual migration 
paths of stars - how long each star spent with its guiding 
centre at each radius. It would predict smaller scale heights 
for populations of stars born in the inner disc. Thus the 
model might predict too high a fraction of high-metallicity 
stars to high altitudes. 

The metallicity distribution at high altitudes depends 
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Figure 7. The model's stellar metallicity distributions at differ- 
ent heights above the plane at Ro. Here we avoided implying a 
specific selection function by using the mass of a specific popula- 
tion to determine its weight in the distribution. The diagrams are 
unsmoothed and the scatter comes from the radial (0.25 kpc) and 
temporal (30Myr) resolution of the model. Upper panel: distri- 
butions of iron abundance. Lower panel: distributions of relative 
oxygen abundances. 
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Figure 8. Upper panel: a scatter plot for a GCS-like measure- 
ment of solar-neighbourhood stars in the ([Fe/H], [O/Fe]) plane. 
Red dots mark positions of stars, while green crosses mark stars 
that are selected to the thick disc via the kinematic selection 
scheme. Lines mark possible criteria to dissect the data with 
a chemical classification scheme in the ([Fe/H] , [O/Fe]) plane. 
Lower panel: the locations in the ([Fe/H], [cn/Fe]) plane of stars 
with spectroscopically determined chemical compositions from 
Bensby (2005), Gilli (2006) and Reddy (2006). 



on the weakly constrained early evolution of the disc and on 
details of mixing, so comparisons with observational data 
would provide valuable constraints on these less secure as- 
pects of the model. Unfortunately, such comparisons are not 
feasible at present . In particul a r, we cannot compare with 
the SDSS data of llvezic et al . (2008) because their metal- 
licity determination breaks down above [Fe/H] ~ — 0.5dex. 
The model however has a constant mean rotational velocity 
in the metallicit y range probed by the SDSS survey, in line 
with the data of llvezic et al ] (|2008h . 

The lower panel of Fig. [71 shows the [O /Fe] distribu- 
tions at different heights. It reveals the bimodal structure 
that motivates the division of the disc into two. The exact 
shape of the two peaks as well as the number of stars in 
between depend on assumptions about gas enrichment and 
the behaviour of SNIa, but the bimodality of the distribu- 
tion is a fundamental prediction of the model, as was shown 
in the appendix of SB09. The increasing bias to high ages as 
| z | increases is reflected in the growing strength of the high 
[O /Fe] peak relative to the low [O /Fe] peak associated with 
the thin disc. 



2.3 The disc divided 

There are principally two strategies by which the disc has 
classically been dissected: by kinematics and by chemistry. 
We caution that different selection procedures do yield in- 
trinsically different samples and that in general these are 
not equivalent. We shall see that these selection differences, 
which account for the spread by almost an order of magni- 
tude in estimates of the relative local densities of the thick 
and thin discs, are readily understood in the context of our 
model. In each scheme criteria are set that define both thin 
and thick disc components, while stars that meet neither 
criterion are here assigned to an "intermediate population" . 
We turn first to chemical selection and then in the light of 
this assess the quality and effects of kinematical criteria. 

The dots and crosses in the upper panel of Fig. [8] show 
a realisation of a GCS-like sample of stars in the model. The 
ridge of the thin disc is evident, as is a ridge of metal-poor 
thick-disc stars at [Fe/H] < -0.65 and [O/Fe] - 0.6. We con- 
sider the thin disc to consist of all stars that lie within the 
black lines around this ridge. Less clear is the extent of the 
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thick disc at [Fe/H] > — 0.6. Guidance is provided by plot- 
ting in green the locations of those stars in the realis ation 
that s atisfy the kinematic selection criteria of iBensbv et al.l 
(2003) to belong to the thick disc. A few of these stars lie 
in the region reserved for the thin disc; this phenomenon 
illustrates the inability of any kinematic selection criteria to 
separate cleanly the thin and thick discs - see 33 below. In 
light of the distribution of green crosses in Fig. [8] we define 
the thick disc to consist of all stars that lie either above the 
horizontal line at [O/Fe] = 0.56 or to the right of the sloping 
line, which has the equation 



[O/Fe] = 0.56 - 0.55([Fe/H] + 0.8) . 



(1) 



The lower panel of Fig. [8] shows the chemical com- 
positions of stars in three large observatio nal programs. 
These studies used different selection criteria - iBensbv et al.l 
_<l2005h kinematically selected for thick-disc stars, while 
iGilli et al.1 (|20Q6h studied stars with planets, so their stars 
are all metal-rich. 

The top and centre panels of Fig. [9] show the distribu- 
tions of rotation velocity (top) and age (centre) within the 
thin disc (green), the thick disc (red) and the intermediate 
population (blue) when the local stellar population is di- 
vided in this way. In the top panel the thick disc stands out 
for the extent to which its V- distribution extends to low V. 
However, its peak lags circular rotation by only ~ lOkms -1 
because it has a significant extension to V > 0. On ac- 
count of its long tail, the average asymmetric drift of the 
thick disc is ~ 22.5 km s -1 , which is slower than that of 
the thin disc by ~ 18kms -1 . The intermediate population 
is much more symmetrically distributed in V and, like the 
V- distribution of the thin disc, peaks near V = with an 
average drift of ~ lOkms -1 . Note that these velocities are 
relative to the local standard of rest (LSR), rather than the 
Su n, which is rotating fast er than the LSR by ~ 5kms _1 
(cf. iDehnen &; Binnevl fl998). Hence relative to the Sun, the 
asy mmetric drift of the thin disc is ~ lOkms -1 . 

IHavwoodl (|2008h showed that the population with mod- 
erate a enhancement is a superposition of stars that either 
combine lower metallicity (at fixed [a/Fe]) with fast rota- 
tion, or higher metallicity with lower rotation. The former 
sub-group bear a clear outer-disc signature, while the lat- 
ter sub-group one associates with the thick disc by virtue 
of their slow rotation. Higher metallicities at given [a/Fe] 
point to a natter trajectory in the ([a/Fe], [Fe/H]) plane of 
the relevant ISM, i.e. to faster metal-enrichment before SNIa 
started to explode. Such fast enrichment is to be e xpected in 
the de nse inner disc. Thus the structure found by IHavwoodl 
(2008|) is an inevitable consequence of chemical evolution in 
the presence of radial mixing. 

The middle panel of Fig. [9] shows that essentially all 
thick-disc stars are older than 6 Gyr. Most stars of the inter- 
mediate population are also this old, but whereas the modal 
age of the thick disc exceeds 10 Gyr, no stars in the inter- 
mediate population are older than 10 Gyr. The sharp rise in 
the number of thick-disc stars at ~ 10 Gyr, just where the 
number of stars in the intermediate population plummets, 
is very striking. The purple curve in the central panel of 
Fig. [10] clarifies the cause of this feature by showing the age 
distribution of stars that have [O/Fe] > 0.56 (the horizontal 
boundary in Fig. [8} and [Fe/H] < -0.8. We see that all these 
metal-poor, highly a-enhanced stars are older than 10 Gyr, 
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Figure 9. The top two panels show the distribution of stars 
from the three chemically selected populations in V velocity and 
age. The bottom diagram shows the distribution of stars by the 
velocity-dispersion parameter of the cohort to which they belong 
- see the description in text. The populations are the thin disc 
(green), the intermediate population (blue) and the thick disc 
(red). The curve showing the age distribution of the thin disc has 
been scaled down by a factor of 10 relative to the curves for the 
other two components. In the other panels each curve is sepa- 
rately normalised to unity. The steps in the age distribution are 
artifacts arising from the model's radial resolution (0.25 kpc); a 
step is produced as an individual ring passes over the selection 
criterion. 



so the significance of a 10 Gyr age is that older stars formed 
before SNIa started to enrich the ISM with iron. The purple 
curve in the top panel of Fig. [10] shows that the modal ro- 
tation velocity of these high- a stars is not far from circular. 
That is, the metal-poor thick disc has a smaller asymmetric 
drift than the portion of the thick disc that lies to the right 
of the division line in Fig. [8] (red curves in Fig. llO[) , which we 
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henceforth refer to as the metal-rich thick disc because all 
its stars have [Fe/H] > —0.8. This metal rich portion has an 
average asymmetric drift of ~ 35kms _1 , which is 30kms -1 
below the mean rotation velocity of the thin disc. 

The green curve in the central panel of Fig. [9] shows 
that all thin-disc stars are younger than 7 Gyr and the rate 
of their formation appears to rise rapidly towards a peak 
at ~ 1.5 Gyr. In reality the SFR in the disc was monotoni- 
cally declining throughout this period, so this apparent rise 
is entirely a selection effect. Several factors contribute to the 
detailed shape of the thin-disc age distribution in Fig. [9] in- 
cluding the restriction of the sample to a volume near the 
Galactic plane (which disfavours old stars) and the bright- 
ening of stars as they begin to turn off the main sequence 
(which accounts for the peak at ~ 1.5 Gyr). 

The bottom panels of Figs. [9] and [10] show decompo- 
sitions of each population into isothermal cohorts. A de- 
composition is possible because in the model each cohort of 
stars (stars formed at a given time and place) has a steadily 
increasing velocity dispersion. The plotted decompositions 
show the distribution of the current velocity dispersions for 
the cohorts that make up each population, weighted by the 
fractional contribution of the cohort to the population. 

In Fig. [9] the isothermal decompositions of the thin disc 
(green) and intermediate population (blue) are similar ex- 
cept that the distribution for the intermediate population is 
shifted to the right by ~ 20kms -1 . The isothermal decom- 
position of the thick disc is bimodal. The purple curve in the 
bottom panel of Fig. [TO] shows that the peak of this curve 
around 40kms -1 is associated with the metal-poor thick 
disc. So within our model the bulk of the metal-poor thick 
disc has smaller velocity dispersions than are found in the 
metal-rich thick disc. This chimes with the higher character- 
istic V velocities of the metal-poor thick disc in indicating 
that it is cooler and faster rotating than the metal-rich thick 
disc. The tail to high dispersions in the decomposition of the 
thick disc is contributed by a small number of very old stars 
that were formed at small radii, where the velocity disper- 
sion is currently large. 

By fitting the model's vertical density profile with the 
sum of two exponentials in \z\, SB09 concluded that in the 
model a fraction /thick ~ 0.13 of solar- neighbourhood stars 
belong to the thick disc; it followed that of order one third of 
the entire disc mass is contributed by the thick disc. These 
numbers were in good agreement w i th th e conclusions that 
ljuric et all (|20Q8h and llvezic et al.1 (|20Q8h drew from SDSS 
counts of stars > 1 kpc from the plane. Using the present 
chemical decomposition into thin and thick discs, we find 
/thick ~ 0.14. In principle this number does not have to 
agree with the value obtained from the density profile. It 
does agree well because at \z\ > 1 kpc the disc is dominated 
by stars that have thick-disc chemistry (Fig. 0. 

Fig. HH shows the vertical density profiles of the thin 
disc (green), thick disc (red) and the entire disc (blue). Fits 
of exponentials to the density profiles yield scale heights of 
268 pc for the thin disc and and 822 pc for the thick disc. 
The two components of the latter have scale heights 690 pc 
for the metal-poor thick disc and 890 pc for the metal-rich 
component. All components show more or less exponential 
profiles. The metal-poor thick disc has the strongest devi- 
ations from an exponential due to its being a mix of very 
old stars from all over the disc with radically different in- 
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Figure 10. As Fig. [9]but with the thick disc split into its metal- 
weak (purple) and -rich (red) parts, the latter being defined to 
comprise [Fe/H] > -0.8. 



trinsic velocity dispersions. When a sum of exponentials is 
fitted to the measured vertical profile of the Galactic disc, 
good fits can be obtained with quite a wide range of scale 
heights on account of a correlation between the scale heights 
of the two discs and their normalisations. The fits above to 
our individual discs are within the range of observationally 
acceptable scale heights (e.g., ljuric et al.ll2008h , consistent 
with the thick disc identified by Juric et al. being close to 
what we have identified in the model using totally different 
criteria. 



2.4 Inside-out formation? 

For simplicity the SB09 model does not accelerate the for- 
mation of the disc at small radii relative to large radii, as is 
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Figure 11. The vertical density profiles of the total disc (blue), 
thin disc (green), thick disc (red), the metal poor thick disc com- 
ponent (purple) and the metal rich thick disc (orange). The thick 
disc partition at z=0 is 14%, scale heights of the single com- 
ponents are (measured between z = 400 pc and z = 3000 pc): 
hthin = 270 pc, h thick = 820 pc, h poorthick = 690 pc and 
hrichthick = 890 pc. 



required by the popular "inside-out" model of disc forma- 
tion. If the model were adjusted to include inside-out for- 
mation, the main change would be to the metal-poor thin 
disc, which would lose parts of its high-1/ wing. Hence the 
inside-out scenario could be put in doubt by demonstrating 
that the V distribution of the high-a stars extends signifi- 
cantly to V > 0, and that many high-y stars have ages in 
excess of 10 Gyr. Further inside-out formation could give rise 
to some alpha enhanced, relatively metal-poor stars younger 
than rsj 10 Gyr by the later onset of star formation in outer 
rings. Neither the thin disc nor the metal-rich thick disc 
would be strongly affected by the introduction of inside-out 
formation. 



3 KINEMATIC DIVISION OF THE DISC 

Because it is much easier to measure the velocity of a star 
than to determine its chemical composition (particularly its 
a-enhancement), nearly all analyses select stars kinemati- 
cally. Our model provides an arena in which we can exam- 
ine the extent to which kinematically selected samples of 
each component will be contaminated with stars from other 
components. 

Samples of local stars such as those of IVenn et all (2004) 
and lBensbv et al.l (2005) are kinematically divided into thin 
and thick-disc stars with the aid of model dist ribution func- 
tions for each component: as described in iBensbv et al.l 
(2003), each star is assigned to the component whose DF 
is largest at t he star's veloc i ty. Bo th DFs are of the type 
introduced bv ISchwarzschildl (|1907l ). namely 



f(U,V,W) = kfiexp 



U 2 (V-Ksym) 2 



2*1 



2al 



W 

^ 2 w 



(2) 



where all components are with respect to the Local Standard 
of Rest, k = (27r) _3//2 (cr {y cry <7v^) _1 is the standard normal- 
isation constant, fi is the relative weight of the population. 
The dispersions a% assumed for the thick disc are larger than 



to be ~ 30kms~ larger for the thick disc than the thin, 
stars with lagging rotation velocities and therefore guiding 
centres at R < Ro also tend to be assigned to the thick disc. 
This effect is reinforced by the fact that the dispersions must 
increase inwards, so stars with guiding centers well inside Ro 
are also likely to be high- velocit y stars . Con sequently the 
"thick disc" stars in IVenn et all (|2004f ) and IBensbv et al.l 
(2005) tend to belong to the in ner disc. In the con text of 
our model this fact explains why IBensbv et al.1 (|2005t ) found 
a long tail of "thick disc" stars that have higher [Fe/H] at a 
given [a/Fe] than the "thin-disc" stars. 

We examine the effectiveness of kinematic selection in 
two ways. First, in each panel of Fig. [12] we plot the distri- 
bution in a Toomre diagram of each of the components that 
we have identified chemically. Subsequently, in Figs [T4HT61 
we examine the distributions in the ([a/Fe], [Fe/H]) plane 
of model stars that have been kinematically identified as 
belonging to the thin or thick disc. 

Fig. 1121 shows the Toomre diagrams for the thin, thick 
and intermediate components, respectively. In Fig. [12] all 
densities are separately normalised to unity, while Fig. [T3l 
shows the density ratios of components in the Toomre dia- 
gram. The extensive overlap of the chemically-selected pop- 
ulations in the Toomre diagram is striking, but a natural 
consequence of the approximately Gaussian nature of the 
distribution functions of each component, which implies that 
the density of thick-disc stars peaks at velocities close to 
the LSR, which is where the thin disc is dominant. Con- 
sequently, no kinematic selection of stars from a particular 
chemical component can be very clean. This point is under- 
lined _bv_th£wlrite curves in Figs. [T2land[T3l which are such 
that Bensbv et all (j2003T ) classified stars witrQ W = 0.5517 
as thick-disc if they lay outside the outermost white curve 
and thin-disc if they lay inside the innermost curve. The top 
panel in Fig. P"2l shows that this criterion does exclude most 
thin-disc stars from a thick-disc sample. However, the up- 
per two panels of Fig.[13]imply substantial contamination of 
the thick disc: in these panels red indicates a region where 
most stars are not thick disc stars, yet at lower right red 
extends significantly beyond the outermost white curve in 
both panels. From Fig.[l2]it is evident that a slightly cleaner 
kinematic separation could be obtained if a non-Gaussian 
distribution function were used in place of J2} , but the main 
problem with kinematic selection is the extensive overlap of 
the components in velocity space. 

From the centre panel of Fig. [12] we see that a large 
fraction of the thick disc is also excluded from a kinemat- 
ically selected sample of thick-disc stars, and many of the 
excluded stars will be assigned to the thin disc because they 
lie within the region reserved for the thin disc. 

Fig. [14] shows the probabilities used by IBensbv et al.l 
(2003) for a star to be assigned to the thin (upper panel) 
and thick (lower panel) discs. The probability of being as- 



1 For general values of W/U the Bensb v et al . (2003) kinematic 
selection criterion, which is three-dimensional, cannot be plotted 
in a Toomre diagram because the latter is two-dimensional. Hence 
we choose the approximate ratio of the dispersion components for 
our graphs. 
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Figure 12. From top to bottom Toomre diagrams for the chem- 
ically selected thin-disc, thick-disc and intermediate populations 
within the model's solar-neighbourhood. Colour encodes the den- 
sity of stars and ranges over 5 orders of magnitude. Shown in 
white are curves on which the probability of star belonging to 
the thick disc by the criteria of Bensby et al. (2003) is constant 
given W = 0.55C7; on these curves from inwards out the thick-disc 
probability is 0.1,1.0 and 10.0 times the thin-probability. Inside 
the inner curve stars were deemed to belong to the thin disc, and 
outside the outer curve they were assigned to the thick disc. 
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V/kms- 1 

Figure 13. From top to bottom the ratio of thick by thin disc, 
thick disc by intermediate disc, and intermediate disc by thin disc 
stars at each point in the Toomre diagram when metal poor and 
metal rich thick disc are combined. Colours show the log of the 
ratio with values in the range (—2.5, 2.5). White contours are the 
same as in Fig. 1121 



signed to the thin disc is large in a broad swath that runs 
from [O/Fe] = 0.6 and [Fe/H] = -1.2 down to the lower 
edge of the populated region of the diagram, and then on 
to solar [Fe/H] and above. Thus the kinematically selected 
thin disc includes high-a stars in contrast to our chemically 
selected thin disc. Kinematic selection does not confine the 
thin disc to stars near the ridge- line of the chemical thin disc 
because the low velocity dispersions and high rotation ve- 
locities characteristic of large radii cause most stars formed 
at large radii to be kinematically assigned to the thin disc. 
Stars in the more metal-rich flank of the chemical thin-disc 
ridge tend to be assigned partly to the intermediate popu- 



lation, or even (for the highest metallicities/innermost rings 
of origin) to the thick disc. 

In Fig.[l4]the probability of a star being assigned to the 
thick disc is high along the sloping upper edge of the popu- 
lated region of the ( [Fe/H], [a/Fe]) plane. Consequently, the 
iBensbv et al kinematic criterion for being a member 

of the thick disc does pick stars that belong to the thick disc 
by our chemical definition. However, the density of stars ac- 
tually assigned to the thick disc by the kinematic criterion, 
which is shown in Fig. [15] extends below the sloping dashed 
line in Fig. [8] because the density of stars assigned to the 
thick-disc is the product of the assignment probability plot- 
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Figure 14. Selection probabilities using the kinematic selection 
function of equation ([2} for the thin (upper panel) and thick disc 
(lower panel). Blue contours give lines of same selection probabil- 
ity for a star at a certain chemical composition with levels running 
from 0.01 to 0.91 with a 0.1 spacing for the thin disc and from 
0.01 to 0.61 with a 0.05 spacing for the thick disc. Colours encode 
the selection probability and the green contours show lines of the 
model's entire disc population density at a 0.5 dex spacing. 

ted in Fig. [TH and the density of stars in the ([Fe/H], [a/Fe]) 
plane, which declines steeply as the edge of the populated 
region is approached. In fact the ridge of kinematically- 
selected thick-disc stars leaves the upper edge of the popu- 
lated region at the a- enhancement turnoff and then runs 
downwards parallel to the thin-disc r idge at an offset o f 
~ 0. 3 dex in [Fe/H], ju st as reported bv lBensbv et al.l ([2003) 
and I Venn et al.l (J200J). Moreover, the zone of thick-disc stars 
merges with the thin-disc popu lation at [Fe/H] ~ , which 
was one of the main findings of Bens bv et al. (2003). 

The upper panel of Fig. [16] shows the age distribu- 
tions of the kinematically-selected components, and should 
be compared with the middle panel of Fig. [9] When kinemat- 
ically selected, the thin disc has a long tail of very old stars. 
Conversely, the age distributions of the thick disc and espe- 
cially the intermediate population extend to much younger 
ages when these components are kinematically selected. It is 
inevitable that a kinematically selected thin disc will contain 
old stars that properly belong to either the thick disc or the 
intermediate population because stars with small velocities 
relative to the LSR must be assigned to the thin disc, yet 




Figure 15. Green lines show the density contours in the 
([Fe/H], [O/Fe]) plane of the entire disc population with a 0.3 dex 
spacing. Colours and blue contours show the absolute density of 
stars selected kinematically according to equation (0 to the thick 
disc with 0.3 dex contour spacing. 



any plausible distribution function for the thick disc will be 
significantly non-zero at such velocities. The assignment of 
young stars to the intermediate population reflects the red 
colour in the lower right corner of the centre panel of Fig. [13] 
that was discussed above. 

The lower p anel of Fig . 1161 shows histograms of the ages 
of GCS stars in iHavwoodl (|2008T ) when stars are assigned 
to com ponents using the kinematic criteria of iBensbv et al.l 
(200 3) - the corre s pondi ng histograms of the ages given 
by iHolmberg et all (2007) is very similar. Clearly the his- 
tograms are badly distorted by errors in the ages, which 
scatter stars to unrealistically large ages, so the horizontal 
scale of the lower panel is nearly twice that of the upper 
panel. None the less, the lower panel seems to be as consis- 
tent with the upper panel as the large errors permit. 

The numbers of model stars in the solar neighbourhood 
that are kinematically assigned to the three components 
analysed in the upper panel of Fig.[16]is thin : intermediate : 
thick = 1 : 0.099 : 0.0239. The same ratios for the obser- 
vational sample analysed in the lower panel of Fig. [16] are 
1 : 0.085 : 0.029 in satisfactory agreement with the model's 
prediction, but the agreement is actually better than this 
comparison suggests, when one accounts for the difference 
between the selection functions used to select the observed 
stars in the lower panel of Fig. 1161 If we use the GCS sample 
without binaries, which our selection function was designed 
for, the ratios are changed to 1 : 0.095 : 0.029. When we 
further remove likely halo stars, the observational thick disc 
fraction shrinks to ~ 0.025. Indeed, the fraction of the lo- 
cal column of thick-disc stars that resides near the Sun is 
sensitive to the distribution of W velocities. The latter is 
not tightly constrained because one of the least satisfactory 
aspects of the model is the absence of dynamical coupling 
between horizontal and vertical motions, which obliges one 
to make an arbitrary assumption about the variation with 
random velocity in the shape of the velocity ellipsoid. It 
is worth noting that the model probably has more metal- 
rich stars high above the Sun than it should have as a con- 
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Figure 16. Upper panel: model age distributions for the three 
kinematically selected components in the model. As in Figures [9] 
and 1101 the curve for thin disc has been lowered by a factor of 50 
and that for the intermediate population by a factor 5 relative 
to the curve for the thick disc. Lower panel: the distribution of 
measured ages of GCS stars given by Haywood (2008) with stars 
kinematically assigned to components using the Bensby (2003) 
criteria. Stars with age younger than 0.5 Gyr are not taken into 
consideration due to high errors. When the ages published by 
Holmberg et al. (2007) are used, the lower panel does not change 
significantly. 



sequence of our assumption that high-velocity stars are as 
susceptible to churning as low-velocity stars. 



4 CONCLUSIONS 

The thick disc is the Proteus of Galactic physics: depending 
on which questions you ask it changes its shape. Although 
it has been identified both by its extended vertical density 
profile and its distinct kinematics, it is most usefully char- 
acterised chemically, not least because chemical composi- 
tion is a permanent feature of a star, whereas distance from 
the plane and peculiar velocity are ever-changing properties. 
Moreover, chemical composition is intimately connected to 
the time and place of the star's birth. 

The determination of the chemical composition of large 
numbers of old main-sequence stars is feasible only for sam- 
ples of nearby stars. Unfortunately, the nearby stars consti- 
tute a strongly biased sample of the whole Galactic disc. It 



is absolutely essential to interpret the statistics of the solar 
neighbourhood in the context of these biases. We have used 
our model Galaxy to explore these biases, and in partic- 
ular the relationship between the components one obtains 
by assigning stars to them on the basis of their kinemat- 
ics or their chemistry. A very straightforward conclusion is 
that kinematic selection inevitably mis-allocates many stars, 
both adding old stars to the thin disc and young stars to the 
thick disc. 

We have shown that our model provides a consis- 
tent interpretation of observations of the solar neighbour- 
hood in which components are identified as regions of the 
([Fe/H], [a/Fe]) plane. Thin-disc stars lie in a narrow ridge 
of high density between [Fe/H] - -0.65 and [Fe/H] - 0.15 
that forms part of the lower edge of the populated part of 
the ([Fe/H], [a/Fe]) plane. The metal-rich thick disc occu- 
pies a broader swath of the ([Fe/H], [a/Fe]) plane that runs 
parallel to the downward-sloping ridge of the thin disc and 
~ 0.3 in [O/Fe] higher. The metal-rich thick disc extends in 
[Fe/H] from ~ —0.9 to well above 0, where it merges with 
the thin disc. At its low-metallicity high-a end, the metal- 
rich thick disc touches the metal-poor thick disc, in which 
[O/Fe] ~ 0.63±0.5 and [Fe/H] goes at least down to - -1.4. 
There is an "intermediate population" of stars that in the 
([Fe/H], [a/Fe]) plane lie between the thin and thick discs, 
but the density of such stars in the ([Fe/H], [a/Fe]) plane is 
relatively low. Thus the two discs are well defined structures. 

Thin-disc stars are all younger than 7 Gyr and are on 
fairly circular orbits. Their values of [Fe/H] and rotation 
velocity V are correlated in the sense that higher V im- 
plies lower [Fe/H]. The stars of the metal-rich thick disc are 
nearly all older than 8 Gyr. Most are on significantly non- 
circular orbits with guiding centres inside Ro and a signifi- 
cant number have V < — 100 km s -1 . Specifically, the metal- 
rich thick disc can be considered to be a superposition of 
isothermal components with radial velocity dispersions be- 
tween 50 and 80kms -1 , strongly peaked around 60kms -1 . 
The metal-poor thick disc consists exclusively of stars older 
than 10 Gyr. Its stars have on average more angular mo- 
mentum and smaller velocity dispersions than the stars of 
the metal-rich thick disc. Among the population of strongly 
a-enhanced stars there is (cf. Fig. 2} an extremely strong 
negative correlation between a and V. 

iMelendez et all (|2008h remarked that the thick disc has 
similar properties to the Galactic bulge. This conclusion is 
natural in the context of our model, in which the metal-rich 
thick disc is made up of stars that have migrated to the 
Sun from the inner disc, where rapid early star formation 
enriched the ISM to significant metallicities before SNIa be- 
gan to lower [a/Fe]. It is to be expected that many of the 
stars that formed alongside the thick-disc stars of the solar 
neighbourhood are now bulge stars. 

Perhaps the most uncertain aspect of the modelling is 
our assumption that a star's probability of being "churned" 
to a different angular momentum is independent of the star's 
random velocity. Since Sellwood & Birmeyj (2002) did not in- 
vestigate the dependence of churning probability on random 
velocity, our assumption could be significantly in error, and 
it is not implausible that stars with large random velocities 
have low churning probabilities. In this case the thick disc 
would be less radially mixed than our model predicts. We 
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will shortly investigate the dependence of churning proba- 
bility on random velocity. 

Given our model's success in synthesising studies of 
the Galactic disc into a coherent picture, it is useful as 
it stands regardless of the theoretical considerations that 
motivated its construction. However, it was not made by 
searching an extensive parameter space for a model that 
would fit the studies described here. Rather it was made 
by building a code that combined standard chemical evolu- 
tion modelling with a model of dynamical evolution that 
refle cts the understandi n g of how spiral structure works 
that Sell wood Sz Binnevl (|2002h gained from N-body mod- 
els and analytical dynamics. Its two free parameters were 
determined from the model's fit to the metallicity gradient 
in the ISM a nd to the metallic i ty dis tribu tion of the GCS 
stars given in INordstrom et al.l (|2004T ) and iHolmberg et al.l 
(|2007T ). Hence we conside r the mod e l's ability to repro- 
duce the data sets of Fuhrmann ( 19981). iBensbv et all (l200l 



Venn et al 



Juric et al 



(120041 ) iReddvet all (120061) iHavwood 
(|2008l ) and llvezic et all (|2008h is remark- 
able and suggests that it has a sound physical basis. 

The key respect in which the model goes beyond tradi- 
tional models of chemical evolution is its inclusion of radial 
migration by stars and inward flow by gas. Inward flow of 
gas is important for the model's success because it estab- 
lishes a much steeper metallicity gradient than traditional 
models produce. The radial migration of stars is absolutely 
key, because it structures the thick disc. Moreover it explains 
the significant spread i n [Fe/H] within the local thin disc, 
and the correlation that iHavwoodl (|2008h identified between 
[Fe/H] and V. 

The discovery that the thick disc overlaps the thin disc 
in [Fe/H] presented a challenge to conventional models of 
chemical evolution because it implies that there are thick 
disc stars that have both [a/Fe] and [Fe/H] higher than in 
some thin-disc stars. The lower values of [a/Fe] in the thin- 
disc stars imply earlier times of birth, so how come [Fe/H] 
is lower? The conventional response to this challenge it to 
suppose that some violent event led to a suspension of star 
formation in the disc, and that during this hiatus a mas- 
sive injection of metal - poor gas lowered [Fe/H] in the ISM 
flChiappini et al. I ll997l . 2001). An objection to this scenario 
is that many other galaxies have thick discs w ith similar 
properties to ours (|Yoachim &; Dalcantonl l20Qq ). so thick- 
disc formation should not require special circumstances. 
We do not press this argument but would strongly make 
the point that a model of chemical evolution that includes 
only essential physics and has a single, early, maximum in 
the star- format ion rate and a monotonically rising value 
of [Fe/H] at each radius automatically produces a thick 
disc with just the properties observed locally. In fact an a- 
enhanced thick disc forms because the speed at which [Fe/H] 
rises declines as one moves outwards, so the value attained 
by [Fe/H] when SNIa start to lower [a/Fe] increases inwards. 
Spiral structure and the Galactic bar scatter a-enhanced 
stars formed at small radii onto more eccentric and more 
inclined orbits and even scatters some of them onto orbits 
of sufficiently high angular momentum that they are found 
in the solar neighbourhood. Readers who want to believe 
in a violent origin of the thick disc may continue to do so. 
But they should be aware that the simplest model of the 
chemo- dynamical evolution of the disc that includes all rele- 



vant physics reproduces the data. Hence there is absolutely 
no evidence that the thick disc has a violent origin. 
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